Effects of Electron-Electron Scattering in Wide Ballistic Microcontacts 
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We calculate the corrections to the Sharvin conductance of ballistic multimode microcontacts 
that result from electron-electron scattering in the leads. Using a semiclassical Boltzmann equa- 
tion, we obtain that these corrections are positive and scale with temperature as T^ \n{Ep/T) for 
three-dimensional contacts and as T for two-dimensional ones. These results are relevant to recent 
experiments on 2DEG contacts. 
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It is a textbook knowledge that normal electron- 
electron scattering does not contribute to the resistivity 
of homogeneous metal. Hovifever if the metal contains mi- 
croscopic inhomogeneities, a quantum-mechanical inter- 
ference between electron-electron and electron-impurity 
scattering results in a temperature-dependent correction 
to the conductivity 1]. The electron-electron scattering 
is also known to affect the resistance of narrow channels 
with boundary scattering because it deflects electrons 
moving along the channel axis to the boundaries, where 
they can dissipate their momentum [3]. Here we ad- 
dress the problem of electron-electron scattering in a con- 
fined geometry without any impurity scattering. Namely, 
we consider a ballistic metal contact with a large num- 
ber of transverse quantum channels and smooth bound- 
aries. Though collisions between electrons do not change 
their total momentum, they change their trajectories and 
hence may prevent some of them from passing through 
the contact or help some extra electrons to get through. 
Surprisingly, to the best of our knowledge this problem 
was never analyzed before. 

In the Landauer-Buttiker formalism the dissipation of 
power in a contact with a perfect transmission is due 
to relaxation processes in the leads [Sj]- This relaxation 
brings the injected electrons in equilibrium with those in 
the electrodes. It is implicitly taken into account by as- 
suming that electrons incident on the contact from the 
electrodes have an equilibrium distribution. However 
this implicit relaxation does not take into account the 
back action of injected electrons upon the electrons in 
the leads, and this is precisely the subject of our paper. 

Studies of ballistic microcontacts were pioneered by 
Sharvin [J]. Generally, these contacts were three-dimen- 
sional (3D) and semiclassical, but both experimentalists 
[5| and theorists [6| focused on the electron-phonon scat- 
tering at high voltages and ignored the electron-electron 
interaction, which was less important in this range. 

In 1990s of the last century, it became possible to fab- 
ricate contacts out of two-dimensional (2D) electron gas. 
The size of these contacts was comparable to the Fermi 
wavelength of electrons, and their conductance exhib- 
ited steps of size e'^/h as the width increased and new 



FIG. 1: The distributions of non-interacting electrons to the 
left and to the right of the orifice jQ]. In the left half-space, 
the Fermi sphere has a dip because the electrons injected from 
the contact have lower energies. In the right half-space, the 
Fermi sphere has a bump because the electrons injected from 
the contact have higher energies. 



quantum channels opened [7|. So the main subject of 
studies became the linear conductance as a function of 
this width. The effects of electron-electron interaction 
were considered for such systems mainly in connection 
with the "0.7 anomaly" [8| in terms of a quasi-localized 
state in the limit of few transverse quantum channels. 
There was also a great deal of theoretical work on very 
long and narrow single-channel contacts with electron- 
electron interactions within the concept of Luttinger liq- 
uid [9]. Very recently, Lunde et al. [lO| studied effects 
of electron-electron scattering in a two-channel long and 
narrow quantum contact using a kinetic equation with a 
collision integral. However no calculations in the limit of 
large channel number have been performed so far. 

We use an approach similar to that of Kulik et al. [a| 
for the case of electron-phonon scattering. As a model of 
the contact, we consider two metal half-spaces separated 
by a thin insulating layer with an orifice. We assume the 
reflection from the insulator to be specular, but all our 
results are valid also for diffuse boundary scattering. The 
radius a of the orifice is much larger than the Fermi wave- 
length but much smaller than both elastic and inelastic 
mean free path of electrons. The distribution functions of 



electrons on both sides of the insulator obey Boltzmann 
equation 

dl dl 
dt dr 

where E — — V(p is the electric field, 
electron collision integral is given by 

d'^k f d'^p' 
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where aee is the dimensionless parameter of electron- 
electron scattering, d = 2 or 3 is the dimensionality of 
the system, and Vd is the corresponding Fermi density of 
states. Equation ^ should be solved together with the 
Poisson equation for the electric potential tp. As the dis- 
tance from the orifice increases, the distribution function 
tends to its equilibrium value fo{£p) and ip tends to eV/2 
and —eV/2 in the left and right half-spaces, V being the 
voltage drop across the contact. The current through the 
contact is given by 
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where p± is the component of p normal to the insulator 
and vector p labels points within the orifice area Sq- 

Because of the condition Ep ^ max(ey, T) one may 
treat the electron velocity near the Fermi surface as en- 
ergy independent and set v = vpn, where n is a unit 
vector in the direction of p. It is possible to avoid solving 
the Poisson equation if one replaces p as the argument 
of / by n and the energy variable e = Ep + e(p{r) — Ep- 
With the new variables, the term with electric field drops 
out from Eq. ([1]), and it takes up the form 



^/(n,£,r) 
dt 



or 



(4) 



Physically, Eq. ^ means that in the absence of col- 
lisions, electrons near the Fermi surface just drift while 
retaining their total energy. The boundary conditions for 
this equation are / = /o(e — eV/2) and f ~ fo{£ + eV/2) 
far from the orifice in the left and right half-spaces. 

Equation ^ may be solved by expanding it in powers 
of aee- In the zero approximation, the electrons simply 
propagate along straight lines while retaining their total 
energy s [ll|, so /(n, e, r) depends solely on whether the 
electron trajectory originates from orifice or not. It is 
convenient to use the notion of solid angle ri(r) at which 
the orifice is seen from point r. In terms of this angle, 
the zero-approximation distribution function is 



/l!!k£>n) = 
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for the electrons in left (upper sign) and right (lower sign) 
half-spaces, respectively (see Fig. 1). A substitution of 
these expressions into Eq. ([3|) results in well known ex- 
pressions for the Sharvin conductance 
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i. e. the conductance quantum times the number of 
transverse channels in the orifice. 

To the first approximation in Uee , the correction to the 
distribution function Sf{n,e,p) at the orifice is readily 
obtained by integrating /ee{/'"^} along the trajectory of 
an electron with momentum direction n that comes to 
point p from infinity 
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Here, r is the time of motion along the trajectory. The 
collision integral in Eq. ([7]) is nonzero only if at least one 
of the momenta in Eq. ([2]) falls within fl{r). As we will 
see below, the main contribution to ([7]) comes from points 
r located much farther from the orifice than a. Hence 
51 (r) may be considered as small and the contribution to 
([7]) from scattering processes where more than one mo- 
mentum lies in 51 (r) may be neglected. It is well known 
that due to the restrictions in the momentum space, the 
electron-electron scattering rate exhibits singularities if 
the momenta of colliding electrons are either parallel or 
antiparallel. The latter singularity [l2| is relevant to our 
case because it results in an anomalous behavior of the 
conductance. The most important contribution to Eq. 
([7]) arises from the collisions of electrons incident on the 
orifice with electrons that are injected from the other 
half-space and have nearly opposite momentum. Hence 
the integration over k in Eq. ^ in the left and right 
half-spaces may be limited to k e 51 (r). The electrons 
with momentum k should be considered as injected, and 
the electrons with the rest of momenta p, p', and k', as 
"native" to the corresponding half-space. 

Now we substitute Eqs. ^ and ^ into Eq. ^ and 
change the order of integrations so that the integrations 
over T and p take place first. The expression for the 
inelastic correction to the current assumes the form 

61 = 2eSQaee'^d^ / ^^ / '^^' / ^^i / ^^^ 
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where the prefactor 2 appears before the integral because 
both half-spaces give equal contributions to SI, the dis- 
tribution dependent factor is given by 

Fl{6,6',6i,62) - [1 - his)] [1 - fnie')] hiei) ^£2) 

-hie) fR{6') [I - /l(£i)] [1 - /l(£2)], (9) 



Q is the Heaviside step function, the quantity 

A(£i, £2, Q) = J ^ J d'^k' d{p' + k' - Q) 

xS{ep,-EF-ei)S{ei,>-EF~e2) (10) 

reflects the limitations in the phase space resulting from 
momentum conservation, and 

f„(n,k) ^l-jd^-^pj dTQ[k e n{p-TVFn)]{U) 

is the maximum possible time of interaction between 
electrons flying to the orifice along the direction n and 
nonequilibrium electrons emitted from the orifice in the 
direction k. As we are interested in the linear conduc- 
tance and the distribution factor F vanishes at V^ = 0, 
we have set (/3 = in all other quantities 

Consider first the 3D case. The quantity fv„ is easily 
calculated based on the geometrical considerations. In a 
spherical system of coordinates with the z axis directed 
along — n, it is given by 

, , , 8ma cos cos 9^ + cos (h sin 6 sin 9^ , , 
T„3(n k) = , . „ =^= , (12) 
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where 6 is the zenith angle of k, (p is its azimuthal angle 
measured from the plane normal to the insulator, and 9^ 
is the angle between — n and the normal to the insulator. 
Note that fVn diverges as 6* — > 0, so electrons moving to 
and from the orifice in the directions — n and n would be 
interacting during an infinitely long time. 

In 3D in the same coordinate system, the phase-space 
factor A is given by 
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where the quantity Ssc — (1/2) max(|ep — ek|, |£i — £2!) 
serves as a natural cutoff for the logarithmic divergence 
that arises in the integral over 9 when performing the 
integration over k in Eq. ([5]). 

Upon integration over k and p in Eq. ([5]) , one obtains 
with logarithmic accuracy 
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The logarithm in this equation may be presented as a sum 
ln{EF/Sec) = \n{EF/T) + ln(T/fec)- As characteristic 
values of Sec in Eq. p4|) are of the order Sec ~ T, the 
second logarithm in this equality may be neglected. Then 
the integrals over the energies are easily calculated, and 
one obtains the relative correction to the conductance 
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In the 2D case, r^ is given by 

_ a cos((/i - 0n) 

T„2(n,k) = ——- — , 
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where </> is the angle between — n and k, and 4>n is the 
angle between — n and the normal to the insulator. The 
phase-space factor A now has singularities at both cf) = Q 
and (p = TT, but only the former is essential for us, so we 
can write it in the form 
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where rj = sin^(0/2) + D/4EJ, and D = [{Sp - £k)^ - 
(ei — £2)^1/4. At _D > 0, the integral of Af,„ over 
logarithmically diverges at = 0, so one has to introduce 
a cutoff Tc = Ic/vp in Eq- (flB)) . The cutoff length Ic may 
be due, e.g., to a very weak electron-impurity scattering, 
which makes the electron distribution isotropic far from 
the orifice. 

Performing the angular integration over and </)„ and 
retaining only the most singular terms, one obtains for 
the correction to the current 

SI2 = — ^ ea^m In I — j de de' / dei / de2 

xS{e + e'-ei-e2)e{D)F{e,e', 61,62)/^. (18) 



The distribution-dependent factor F in (jlSp may be lin- 
earized with respect to V. With account taken of the 
delta function in this equation, 

F(£,£',£i,£2) = ^e^/o(£)/o(e')/o(ei)/o(e2). 

Then (J18p is easily evaluated by introducing dimension- 
less variables ^i — e.i/T, and we arrive at an expression 
for the relative correction to the conductance 
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In both 2D and 3D, the corrections to the conduc- 
tance are positive, unlike the cases of electron-phonon 
and electron-impurity scattering. This is due to a prin- 
cipally different nature of electron-electron scattering, 
which does not result in a momentum relaxation of elec- 
trons. The positive contribution to the conductance may 
be considered as a sort of drag of electrons incident on 
the contact by electrons emitted from the other elec- 
trode with a different velocity. As shown in Fig. ^ 
the electron-electron scattering tends to smooth out the 
bump in the angular distribution of electrons in the right 
electrode while conserving their total momentum at a 
given point, so the number of electrons at the opposite 
side of the Fermi surface decreases. In the case of a dent 
in the left electrode this number increases and hence the 
two corrections add together. 




electron scattering considered here. For the electronic 
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FIG. 2: Scattering processes that result in a positive correc- 
tion to the current in 2D. The gray arrow shows the direction 
of net particle current. Extra electrons above the Fermi sur- 
face (black circles) are injected from the contact and reside 
within the solid angle 0(r) (shaded area). As an extra elec- 
tron in state k collides with electron p, they are scattered into 
k' and p'. Momentum conservation requires that p' = — k' 
and p = — k, hence the relaxation of extra electrons decreases 
the number of particles that will hit the orifice. 



Naively, one might expect 5G/Gq to be of the order of 
a/lee, where l~^ ~ vpT'^/Ep is the inverse equilibrium 
electron-electron scattering length. However the actual 
effect is much larger: 6G3 and SG2 contain additional fac- 
tors \n{Ep/T) and Ep/T, respectively. These factors are 
due to a sort of " resonance" between electrons with oppo- 
sitely directed velocities moving to and from the contacts, 
which results from a superposition of the singularities in 
the scattering rate and the effective time of interaction of 
these electrons. The anomalously large correction to the 
conductance is closely related to the logarithmic singu- 
larity in the lifetime of quasiparticles in a 2D electron gas 
[13] but is enhanced by the nonequilibrium distribution 
function with a strong angular dependence. The mag- 
nitude of the effect is increased because the scattering 
intensity falls off away from the orifice so slowly that the 
actual size of the scattering region is much larger than 
the diameter of the orifice. 

Very recently, Renard et al. [14] observed that the 
conductance of wide 2D quantum contacts with a large 
number of transverse channels increased with tempera- 
ture. The increment was proportional to the tempera- 
ture and contact width. The authors attributed this in- 
crease to a quantum-mechanical interference between the 
scattering of electrons from impurities and from Friedel 
oscillations around them 



is r^ ^ 1 and kp is of the order of inverse screening 
length, hence it is reasonable to take a^^. ~ 1 1]. With 
a ^ Af , we get the coefficient of T in ( ^ of the order of 
XQ^^K^^, which roughly agrees with 14| . 

The authors of [IJI also report an anomalous positive 
magnetoresistance in low fields. This effect may be easily 
explained in terms of our results because the magnetic 
field bends trajectories of oppositely moving electrons 
in opposite directions and destroys the "resonance" be- 
tween them. The positive correction to the current from 
electron-electron interaction should be suppressed when 
the classical cyclotron radius Ic = pp/eB becomes of the 
order of the size of the interaction region aEp/T. 
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ment may be due to purely classical effects of electron- 
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